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Abstract. Numerical and physical experiments on two-dimensional {2d) turbulence show 
that the differences of transverse components of velocity field are well described by a gaussian 
statistics and Kolmogorov scaling exponents. In this case the dissipation ffuctuations are 
irrelevant in the hmit of small viscosity. In general, one can assume existence of critical 
space-dimensionality d = dc, at which the energy flux and all odd-order moments of velocity 
difference change sign and the dissipation ffuctuations become dynamically unimportant. At 
d < dc the flow can be described by the "mean-field theory" , leading to the observed gaussian 
statistics and Kolmogorov scaling of transverse velocity differences. It is shown that in the 
vicinity of d = dc the ratio of the relaxation and translation characteristic times decreases 
to zero, thus giving rise to a small parameter of the theory. The expressions for pressure 
and dissipation contributions to the exact equation for the generating function of transverse 
velocity differences are derived in the vicinity of d = dc- The resulting equation describes 
experimental data on two-dimensional turbulence and demonstrate onset of intermittency as 
d—dc > and r/L — > in three-dimensional flows in close agreement with experimental data. 
In addition, some new exact relations between correlation functions of velocity differences 
are derived. It is also predicted that the single-point pdf of transverse velocity difference in 
developing as well as in the large-scale stabilized two-dimensional turbulence is a gaussian. 
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1 Introduction 



The role of the mean-field theories and gaussian limits as starting points for understanding 
of such important physical phenomena as superconductivity, superfluidity, critical point, 
naming just a few, can hardly be overestimated. These theories, usually based on remarkable 
physical intuition and insight, provided mathematical and intellectual basis for investigation 
of much more difficult regimes in terms of deviations from the mean field solutions. The most 
recent example is a theory of anomalous scaling in a model of a passive scalar, advected by 
a random velocity field, which was developed as an expansion in powers of small parameters 
characterizing deviations from the two gaussian limits [l]-[3]. In a typical non-linear system, 
a gaussian limit corresponds to a weak coupling asymptotics and, as a consequence, to a 
"normal" , non- anomalous, scaling, which can often be obtained from a "bare' or linearized 
problem. A good example of this behaviour is a fluid in thermodynamic equilibrium. 
The large- Reynolds- number three-dimensional {3d) strong turbulence is characterized by an 
0(1) energy flux S — i'{diVjY, which in many fiows is 0{vf^g/ L), where L is an integral scale 
of turbulence. In the inertial range, where ka » /c— >oo or r/ L « 1, the observed energy 
spectrum E{k) is close to the one, proposed by Kolmogorov and the probability density 
P{Au) with Am = u{x + r) — u{x), is far from the gaussian. Moreover, the experiments 
revealed that the moments of velocity difference Sn,o = (Am)" oc r^" with the exponents ^„ 
given by very strange ("anomalous") numbers, which cannot be obtained on dimensional 
grounds. This anomalous scaling and the very existence of the energy fiux, resulting in the 
non-zero value of the third-order moment -S'a^o — (Am)^ ?s 0(r), where u • r = ur, implies a 
strongly non-gaussian process and an obvious lack of the mean- field limit. 
The situation may not be so grim, however: all odd-order moments of transverse velocity 
differences in both 2d and So? flows 5'o,2n+i — ('^ (^^ + r) — v{x)y^+^ — with v • r = 0. This 
fact tells us that these components of velocity differences do not participate in the inter-scale 
energy transfer and there is no a'priori any reason for them not to obey gaussian statistics in 
some limiting cases. This is indeed true in two-dimensional turbulence in the inverse cascade 
range where If « r « L and // is a forcing scale. 
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Numerical and physical experiments on the external- force- driven two- dimensional turbu- 
lence showed that the moments of transverse velocity differences and even-order moments 
of the longitudinal ones are very close to the gaussian values and are characterized by the 
Kolmogorov scaling exponents <5'o,2n oc S2n,o oc r«2" with = f [4]- [7]. The odd-order 
moments of longitudinal velocity differences are positive in a 2c? flow, while they all are 
negative in 3d. This observation tells that the most important distinction between two and 
three-dimensional turbulence in in dynamic role of the dissipation contributions: they are 
irrelevant in the two-dimensional inverse cascade range and are crucial (see below) for the 
small-scale dynamics of a three-dimensional flow where the forcing terms can be neglected 
(see below). Thus, we can assume that it is the dissipation fluctuations that are responsible 
for both strong deviations from the gaussian statistics and anomalous scaling in three di- 
mensional flows. This assumption is consistent with the observation of the close-to-gaussian 
probability density of velocity differences in 3d turbulence at the scales r ^ L where the 
integral scale L is defined as the one at which: 5*3^0 (-^) — and the intermittent dissipa- 
tion fluctuations disappear [8]. It is clear that this range (r ~ L) is not characterized by a 
well-defined scaling exponents. 

As will become clear below, the forcing contribution to the equation for the probability 
density involves factor ji ^ 1 — cos{kfr). This means that at the scales r >> 1/kf the 
parameter /i ^ 1 while /i ~ {kfrY when kfV « 1. Thus, the forcing term must be 
important in the inertial range of two-dimensional turbulence (r >> //) and is irrelevant 
in the 3(i-inertial range with the positive energy flux. This change happens at some space 
dimensionality d — dc 8ki which the energy flux changes sign. First calculation of dc ~ 2.05 
was conducted by Prisch et. al. [9] within the framework of a simple closure model. The 
more physically transparent calculation can be performed for the Navier-Stokes equations 
driven by a random force having an algebraically decaying spectrum in the inertial range [10], 
where a one- loop small-scale- elimination procedure gives a correction to the bare viscosity: 

d^ -d-e 2 
dv — V — -tie 

2d{d + 2) 

where Re is a properly deflned Reynolds number corresponding to the eliminated small-scale 
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velocity fluctuations and e is a parameter characterizing the forcing function. In case of 
Kolmogorov turbulence e ~ 4. This relation shows that the role, the small-scales play in 
turbulence dynamics, depends on the space dimensionality d: the correction to viscosity is 
positive when d > dc{e) where it changes sign. Physically, this means that the small-scale 
velocity fluctuations take energy from the large- scales motions (direct energy cascade from 
large-to-small structures) ai d > dc, while at d < dc they excite the large scale motions (spend 
their energy) giving rise to the inverse energy cascade. For e — 4 the critical dimensionality 
dc ~ 2.56. The correct value of dc is not too important: what is crucial for the theory 
presented below is that the critical dimensionality, at which the flux changes its sign, exists. 
It will be shown below that d — dc ^ is a. small parameter of the theory enabling one to 
calculate an expression for the dissipation anomaly in the form, resembling the Kolmogorov 
refined similarity hypothesis. 

Since in 2d the moments 5'o,2n show Kolmogorov scaling, the gaussian statistics of transverse 
velocity differences cannot correspond to the weak coupling limit. This problem was consid- 
ered in the Ref. [11] where, following Polyakov [12], the equation for the generating function 
for the problem of the Navier-Stokcs turbulence was first introduced . An unusual symmetry 
of this equation enabled one to show that the solution was consistent with both Kolmogorov 
scaling and gaussian statistics. In this work a more detailed theory of two-dimensional tur- 
bulence is presented and the generalization to three-dimensional flows is considered. The 
main result of the paper is a model demonstrating how the deviations from the "normal 
scaling" and gaussian statistics appear in 3d when the strength of the dissipation term /3 > 
and the "scale" parameter e — 1 — r/L deviate from zero. 

This paper is organized as follows. In Section 2 the equation for the generating function, 
derived in [11] is introduced. Some new exact relations between velocity structure functions, 
following from this equation, are derived in Section 3. The connection between scaling expo- 
nents of the moments of velocity differences and their amplitudes is established in Section 4. 
The mean-field derivation of the pressure term is given in Section 5 which is used to obtain 
a gaussian pdf in the two-dimensional flow in Section 6. In Sections 7 the small parameter 
of the theory is identifled and used for derivation of the expression for the dissipation con- 
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tributions to the eqaution for the pdf. Section 8 is devoted to solution of the equation and 
demonstration how anomalous scaling and deviations from gaussian statistics emergr from 
the theory. Conclusive remarks are presented in Section 9. 



2 Equation for the Generating Function 

The equations of motion are (density p = 1): 

dtVi + VjdjVi ^ -dip + uW^Vi + fi] diVi^O (1) 
where f is a forcing function responsible for the large-scale kinetic energy production and 



in a statistically steady state the mean pumping rate P — f ■ v. In what follows we will be 
mainly interested in the probability density function of two-point velocity difference U = 
u(x') — u(x) = Au. The generating function is: Z —< exp{X • U) >. The equation for the 
generating function of velocity differences corresponding to (1) is: 

with 

7/ =< A • Af e^-^" > (3) 
= -A- < e^-^"A(Vp) -A- < e^-"(V2p(a;2) - Vip(a;i)) > (4) 

and 

D = u\- < (V^v(x2) - V?v(xi))e^-u > (5) 

The most interesting and surprising feature of (2) is the fact that, unlike in the problem of 
Burgers turbulence [12], the advective contributions are represented there in a closed form. 
This means that the theory, developed below, is free from the troubles related to Galileo 
invariance, haunting all schemes, based on renormalized perturbation expansions in powers of 
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Reynolds number. To completely close the problem the expressions for Ip and D are needed. 
The equations (2)- (3) formulate the turbulence theory in terms of "only" two unknowns Ip 
and D. The Kolmogorov refined similarity hypothesis stating that (Au)^ — (t)£rr where 
is a scale-independent random process and £r is a dissipation rate averaged over a ball of 
radius r around point x, can be a promising starting point to a closure for the dissipation 
term D. This will be done below. The pressure term in (2)-(3) is also of a very specific and 
rather limited nature: all we have to know is the correlation functions < UiUj ■ ■ ■ C/,„AVp >. 
Thus, the definite targets needed for derivation of the closed equation for Z-functions are 
well-defined. 

The generating function can depend only on three variables: rji — r; ri2 = ^ = 

Xcos{9); ri3 = \J \^ - rjh and 

Zt + [dA + ^5.. + jdr,,d,, + ^^^£^0,, - ^-^dl]Z ^If + Ip + D (6) 

where, to simplify notation we set di^a = and v{i) = f(xi). Below we will often use 
= di. The functions Ip, If and D are easily extracted from the above definitions. Let us 
denote Au = U and Av = V. 

In the new variables the generating function can be represented as: 

with the mean dissipation rate £ defined by: = ^£. Any correlation function is 

thus: 

Sn,m =< U^V^ >= d^d^Z{ri2 = ?73 = 0, r) 

3 Relations between moments of velocity difference 

Let us discuss some direct consequences of equations (l)-(6). The Navicr-Stokcs equations 
are invariant under transformation: v — » —v and y —y. That is why: < {dyp{0) — 
dyp{r)){Av)"^ >7^ if m = 2n -|- 1 with n > 1 and is equal to zero if m = 2n. It is also clear 
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from the symmetry that < dxpU'^"' >= 0. It follows from the Navier-Stokes equations that 

< (V^C/)^/^" > = < (V2l/)l^2n 0. 

Multiplying (6) by 773, applying (?|"~^ to the resulting equation gives as 772 = 773 —> 



+ '52n,0 '52n-2,2 = - COs{r/lf)) 52n-3,0 

(7) 

The right side of (7) is O(r^) in a three-dimensional flow (3d) where r/lf « 1 and and 
can be neglected. It is however is 0(1) in two dimensional turbulence in the inverse cascade 
range where r/lf » 1. The dissipation terms do not contribute to this relation. For n = 1 
(7) gives a well-known incompressibility relation [13], [14]: 

9-^2,0 .d-'^c^ _d-l . . 

1 >J2,0 — '-'0,2 [o) 

Multiplying (6) by 773 with subsequent d^&l leads after setting 7^2 = % = as 1/ — > to: 

^ + ^53,.-2^^5,, = (-l)<'ip (9) 

or r r a 

where d = 2; 3 . Applying 9|7;3 to (6) gives, as — > 0: 

Substituting this into (9) yields a well-known Kolmogorov relation: 

S,„.(Auf = (-ir^^Pr (11) 
For 271 = 4 the relation (7) reads: 

dS,nd-l _ 3(d-l) 

1 -^4,0 — '-'2,2 U^j 

or r r 
In two-dimensional turbulence {d — 2) in the inverse cascade range one can neglect the 

dissipation contribution D (see below) and derive: 

+ ^^5i,2n = n{2n - l)PSo,2n-2 - 2n < Vy.{^vf^-' > (13) 
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where Vyv = dyv{x + r) — dyv{x). Another interesting relation, vahd in 2d is obtained from 
(6) readily: 



+ ^^S2,2n = -2n< Vy^AuiAvf--' > +n{2n - l)PS^,2n-2 (14) 

In the direct cascade range, where the forcing contribution is O(r^) —>■ 0, the relation (14) 
for an arbitrary dimensionality d reads: 

dS2,2n ^d-l + 2n 2n + d-l 5o,2n+2 o / -r? a ^ a ^2n-l ^ n rx 

a H S2,2n = -i: — — ; ■ 2n< VyvAu{Av) " > (15) 

dr r 2n + l r 

Now, let us multiply (6) by 773, differentiate once over 772 and three times over 773. This gives: 

882,2 d + d+1 



+ 82,2 = ^—^0,4 - 2Vy,AuAv (16) 

or r Sr 



This relation is correct since uV^vAuAv — uV^u{Av)'^ — 0. 

2d simulations of Bofetta, Celani and Vergassola To achieve a true steady state these 
authors [7] conducted a series of very accurate simulations of the problem (1) with the 

large-scale dissipation term Di = — av in the right side (6). The moments of transverse 
velocity differences, reported in this paper (n > 2), were very close to their gaussian values. 
It is clear that this term introduces —aX^-^ into the right side of (6) which is small in 
the inertial range where the non-linearity is large. One has to be careful though with the 
dangerous interval A^; — > where the linear term is not small. We expect the negative- order 
structure functions with — 1 < ri < strongly depend on the functional shape of the otherwise 
irrelevant large-scale dissipation term. The same can be predicted for various conditional 
expectation values of dynamical variables, like pressure gradients and dissipation terms, for 
the fixed values of AvAu: near the origin where Av and Au are very small, the artificially 
introduced linear contributions to the Navier-Stokes equations dominates, producing large 
and non-universal deviations from the universal functions characterizing inertial range. For 
example, with addition of the linear dissipation, the relation (14) reads: 
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(17) 

modifying the balance (pressure contribution) in the range of small product Am or r/L fa 
1. In the interval where AuAv is not small (r — > 0), the linear terms are small and can be 
neglected. 

The results obtained in this section can also be derived with 

Z = e^^''3^+''2^ (18) 
with the properly defined moments of V and U. 



4 Asymptotic values of exponents in three-dimensional 
flows 

In case of intermittent turbulence -S'^^^ = 74^_„r^'"'". We can see that in the inertial range 
of a three-dimensional fiow {r/lf « 1) the right side of (7) is neghgible and, as a result, 
^2n,o — ^2n-2,2 = ^271 < Substituting this into (7) gives immediately: 

6n = (ci-l)[(2n-l)^%^-l] (19) 
Let us introduce the probability density functions: 

S2n,0 ^U^^ j P{U)U^^ dU (20) 

and 

S2n-2,2 = ^2"-2y2 = J p{U)U^''-\{V\U)dUdV (21) 

where q{V\U) is conditional pdf of V for fixed value oiU. It is clear that q{V\U) — q{—V,U), 
so that all odd-order moments of V are equal to zero.. This expression can be rewritten as: 

S2n-2,2 = = J P{U)U^"-^Q2{U)dU (22) 
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where Q2 is a conditional expectation value of V"^ for a fixed value of U : 



Q^iU) = j V\{V\U) dV (23) 

Comparing the above expressions we observe that the amphtudes A2nfl ~ ^2n-2,2 only if 
Q2 oc U"^. This seems rather improbable in the limit U 00 {n ^ 00). As a result, the 
linear regime ^2n oc n is equally improbable. 

Saturation of exponents ^2n ^00 — const as n ^ 00 is possible on a rather wide class of 
probability densities. For example, 

P{U) ^ -A^^P(U)Q2{U) (24) 
where ^ > is a constant. Then, 

S2nfi = (2n - l)A j P{U)Q2{U)U^''-''dU = (2n - l)AS2n-2,2 (25) 

Substituting this result into the expression for ^2n gives ^2n {d — 1){A — 1) — const. The 
relation (24) defines the large-C/ asymptotics of the pdf P{U) in terms of Q2{U): 

1 A udu 

P{U)<x-^e-^J ^ (26) 

HQ 2(17) — > U^, then, assuming the existence of all moments, it follows from (25) that /3 < 2. 
The "lognormal " pdf P{U) corresponds to Q2{U) oc U'^/2log{U). The expression (24) also 
gives in the limit of large n: 



S2n+1,0 = = 2n>ly2[/2n-l ^ 2nAS2n-l,2 (27) 

5 Pressure contributions 

Due to the symmetries of the Navier-Stokes equations, neither pressure nor dissipation terms 
contributed to the expressions (7) -(12). To proceed further we have to evaluate Ip and D. 
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First of all we see from (12) that ,^4^0 = ^2,2- Let us assume that in the inertial range 
'S'4,0 = A^^r^^-o, 5*0,4 = ^0,4?^^°''' and ^2,2 = ^2,2r^2'^ Then, it is clear from (12) and (16) that 
neglecting the pressure contribution to (16) gives ^4,0 = Co,4- 

d = 2. It will be shown below that in 2d the even-order moments of velocity differences 
are very close to the gaussian ones and all exponents are close to the K41 values ^„ = n/3. 
Then, A^^ = 3A|q and Ao,4 = 3Aq2. It follows from (12) that ^2,2 = |^4,o- Taking into 
account that when d — 2 the amplitudes §^2,0 = ^0,2, we conclude that without the pressure 
contribution the equations ( 12) and (16) are incompatible. 

Following [15] we introduce a conditional expectation value of the pressure gradient difference 
for a fixed value of Au, Av and r: 

< dyp{x + r) - dyp{x)\Au, Av, r>^Yl «m,n(r)(A«)'"(A^;)" (28) 

m,n 

where the functions Km,„(r) ensure proper dimensionality of the corresponding correlation 
functions. The above expression explicitly assumes existence of an expansion of conditional 
expectation value (28). In general, this may not be true due to various singularities such 
as the ones arising in the dissipation contributions (see below). Since the pressure term 
involves only one spacial derivative, the ultra-violet singularity cannot appear. The infra- 
red singularity is not there at least in 2d where the integral scale is time-dependent (see 
below). Keeping only the first two terms of the expansion (28), produces a model for the 
pressure contributions: 

^^zi^^v Av 

< dyp{x + r) — dyp{x)\AyAv >^ —h b j (29) 

r (Pr)3 

Since in an incompressible and homogeneous fiow AuVxu — AvVyv — 0, the coefficients h 
and h are related as: 

-/i5i,2 = 65o,2(Pr)3 (30) 

Limiting the expansion of a conditional expectation value by the first terms resembles Lan- 
dau's theory of critical phenomena, well describing experimental data in a certain range of 
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parameters variation. We will show below that in case of turbulence this approximation 
gives the results which are in agreement with the data. This may be a consequence of the 
fact that Ao,4 ~ ^4,0 ~ ^2,2 — 0{1). 

With — n/3 it follows from (7), (13)-(15) and (29) that when d — 2 and n — > 00: 

( hi) =n{2n-l)PSo,2n-2 + b 



(Pr) 



and 



^MA-Shl ^^^S2^^ ^ - 1)PS, 2n-2 + 

3 r r ' (Pr)3 



In the limit n ^ 00 assuming that Si^2n ~ nSi^2n-2 one derives readily: 

PrnSo,2n-2 ~ ~ p^So,2n+4 (31) 

which is consistent with the gaussian pdf P{Av) as Av — > oo. Thus, the relation (29) implies 
the gaussian tails of the probability density. It is clear that due to the finite energy flux 
and relations (lO)-(ll), two-dimensional turbulence cannot be a gaussian process. All the 
relation (31) can tell us is that the even-order moments with n » 1, described by (31), can 
be close to the gaussian values. It will be shown below that transverse velocity differences, 
not directly involved in the inter-scale energy transfer, can obey gaussian statistics. 
It is clear from (17) that the model (29) for the pressure contributions is wrong when the 
linear dissipation terms are added to the Navier-Stokes equations. In the limit of small 
AuAv the balance is achieved when: 

^ I A A , AuAv , b , ^ 

< VyvlAu, Av, r >Ri -h (a H 2-)Af 

r (Pr)3 

which differs from (29) in the range of small AuAv. 

3d. In the intermittent three-dimensional turbulence ^2n < C2n+i- This produces strong 
restrictions on the structure of the pressure contributions to the equation (6). Let us assume 
that ^2n,o = C2,2n-2 = ^2n-2,2- Then, it is clear from (15) that the first term of expansion (28) 
has all right properties. The relation (15), involving the r-derivatives, is valid for an arbitrary 
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value of n and that is why any additional term of expansion (28) must not only depend on 
a proper power of n but the functions ^^^^(r) must also reflect non-trivial dimensionalities 
caused by the anomalous scaling exponents It cannot be rigorously ruled out, though 
this possibility seems quite bizzarre. In what follows we will adopt the pressure model (29) 
in the three-dimensional case also. 

6 Two-dimensional turbulence 

Now we will be interested in the case of the two-dimensional turbulence in the inverse cascade 
range. If a two-dimensional (2d) fluid is stirred by a random (or non-random) forcing, acting 
on a scale = 1/kf, the produced energy is spent on creation of the large-scale (/ > //) flow 
which cannot be dissipated in the limit of large Reynolds number as u ^ 0. This is a direct 
and most important consequence of an additional, enstrophy conservation, law, characteristic 
of two dimensional hydrodynamics [16]. As a result, the dissipation terms are irrelevant in 
the inverse cascade range and we set D — in (6) and hope that in two dimensions the 
situation is greatly simplifled. This hope is supported by recent numerical and physical 
experiments [4]- [7] showing that as long as the integral scale (x.ti is much smaller than 
the size of the system, the velocity held at the scales Lj >> Z >> is a stationary close-to- 
gaussian process characterized by the structure functions with the Kolmogorov exponents 
— iT'/S. In a recent paper Bofetta, Celani and Vegrassola [7] reported the results of very 
accurate numerical simulations of two-dimensional turbulence generated by a random force. 
No deviations from gaussian statistics of transverse velocity differences as well as from the 
Kolmogorov scaling ^„ = n/3 were detected, and is equal to zero if m = 2n. 
The pressure gradient dyp — dydidjd~'^ AviAvj and the difficulty in calculation of Ip is in the 
integral over the entire space defined by the inverse Laplacian d~'^. The huge simphfication, 
valid in 2d, comes from the fact that all contributions to the left side of equation (6) as 
well as // are independent on time. This means that the integrals involved in the pressure 
terms cannot be infra-red divergent since in a two-dimensional fiow L — L{t) oct^. We also 
have that ip — > a^ip when U, V ^ aU ; aV. Based on this and taking into account that 
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< (At')^"'+^(M^ + '^y + UyVx) >= 0, we, in the limit r}2 — > 0, adopt a low-order model (29) 
giving: 

In two dimensions the relation (32) combined with (6) in the limit 772 ^ 0, solves the problem 
2d turbulence. 

Substituting (32), (29) into (6) and, based on (9)-(ll), seeking a solution as 772 —> as 



1 111 Tj^ 

Z{m,V3,r) ^ Z3{r)3,r)ip{r)2r^,r)2) « Z2{r)3,r)exp{-A2,o{r)2Pr^f){l+-Ai^2vh2{Pr)+^+...) 

(33) 

where (^1,2 = 1/2), gives: 

[dr + - + ^Vsds] ^^vIZs = 2Pr,IZs + ^d,Z, (34) 
r r 2 r3 

Setting Z3 = Z^irisrs) = Zj,{X) and /i = 4/3 one derives using, the relation (30) (6 = -3^): 



2^0,2-^^ -^3 — QxZ-i 

corresponding to a Gaussian solution with the correct width /lo,2- This fact serves as a con- 
sistency check that the gaussian is a solution for the PDF of transverse velocity differences. 
The equation (34) defines a probability density function corresponding to the finite moments 
'S'2m,2n('^) Only whcu /i = 4/3. This situation resembles Polyakov's theory of Burgers turbu- 
lence [12] reduced to an eigenvalue problem with a single eigenvalue corresponding to the 
pdf which is positive in the entire interval. 

Having these exact results, and keeping in mind (33) one can integrate the equation over 772 
from — ioo to to obtain : 

— h3(l — At — 0) — - — = r^s^s (35) 

dr ^ ' r dr)3 (Pr)5 

valid as long as 
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0^2,0 

This constraint is an artifact of an approximate relation (33). As will be shown below 
(35) gives an exact gaussian solution and thus, is valid beyond above interval. This result is 
obtained choosing the integration function ^(7^3, r) to compensate the 0{Z/r) term violating 
the normalizabihty constraint Z(0, 0,r) — 1. Solution to (35) is: 

Z3^exp{^r]l{Pr)l) (36) 

with the parameter 7 = oc Aq^2 defining the width of the gaussian. 

The first-order differential equation (35) for the generating function differences imphes the 

underlying linear Langevin dynamics of transverse velocity differences. It is important that 

this equation in non-local in physical space but local in the Fourier one. The effective forcing, 

corresponding to the right-side of (35), is a non-local and solution- dependent. 

To evaluate the single-point probability density, corresponding to velocity differences in the 

limit of large displacements r , we can notice that the energy flux is not equal to zero only 

at r << L. At the distances r > L{t) the zero value of the energy flux and symmetrization 

of the probability density {P{Au, L) = P(— Am, L) can be achieved only when the pressure 

contribution to the equation (6) compensates the advective terms. As a result, since D — 0, 

we have: 

Zt - 2Pr]lZ 

Seeking a solution ets Z — Z{rj3y/i) = Z{r]sVrms{t)) gives a gaussian result: 

Similar outcome is obtained for the case investigated in [7]. When turbulence is stabilized 
at the large scales by an artificially introduced friction, the resulting equation is: 

rjadaZ = 2pirjlZ 
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also leading to the gaussian pdf. To conclude this section we would like to discuss the 
physical meaning of the integral scale L. The integral scale of turbulence is a scale at which 
the flux decreases to zero [18] and at which 5'3,o(-^) = 0. 



7 Small Parameter in Turbulence Theory in Three 
Dimensions 

Three-dimensional turbulence is a notoriously difficult problem due to absence of the small 
parameter. This can be illustrated on an example of a coarse-graining procedure, which was 
extremely successful in engineering turbulence simulations. Consider the wave-number k in 
the inertial range. Let us denote v^(k) the Fourier components of the velocity field with all 
modes v^(q) = when q > k. The coarse-grained field in the physical space is defined then: 

v^(x,t) = /' ^ Y<{k)e^^''d^k 

r 

The equation of motion for vf{k) in the Fourier space resembles the Navier-Stokes equation 
with effective viscosity 

with dc ~ 2.56 plus high order non-linearities. The parameter 3.2 in the above relation, 
evaluated at d = 3, is in fact a weak function of space dimensionality and S = P = 0(1). 
This expression is derived assuming a close-to-gaussian statistics of the small-scale turbulence 
with 

Viik,uj)vAk' ,uj') oc — : ——-T5{k + k')5{uj + uj') 

—luj -\- ^{kjk^ 

which is accurate a,t d — dc ^ [10] (see below). 

In the physical space the effective viscosity of the coarse-grained field: 

Ur^{d- dc)^N{Sr^)^ ^ v'^Tr 
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defines the relaxation time which is a characteristic time of interaction of the field Vr with 

the eliminated modes acting on the scales I < r. 

The difficulty of the theory is in the higher non-linearities 



and the dimensinless expansion parameter 



Tr _ TrVr 
v' r'^r^v 



Br T 

is nothing but the ratio of translational time Or, characterizing the tendency of the "large- 
scale" longitudinal velocity fluctuation at the scale r to form a small-scale structure ( "shock" ) 
in the absence of pressure, to the relaxation time strongly influenced by the pressure gradients 
contributions. In both 2d and 3d these times are of the same order and that is why trancation 
of the expansion is a very difficult problem. 

This is not so in the vicinity oi d — dc- Let us assume that the theory can be analytically 
continued to the non-integer dimensions [9]. Then, since the energy flux is 0{1) we have: 

£ - 1- < A^A.)^ >^^< {Aur . < i^r >= 0(1) 



dr dr drj 



This means that 



Vr,rms ~ {d - dc) ^{Et)^ 

the dissipation wave- number ~ {d — d^^{-^)^ and, as a consequence, 

J ~ (c?-4)^ 

Or 

which will serve as a small parameter of the theory when d — dc^ 0+. 

These relations tell us that the turbulent intensity grows to inflnity with d — dc ^ where 

the energy flux changes its sign. The time needed to reach the steady stae is estimated 

easily: 
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1 12 

T {d — dc)~^S~^r3 

after which a close-to-Kolmogorov spectrum is expected both above and below dc- Thus, at 

d = dc the the flow is unsteady. 

The above results enable one to derive a plausable estimate for the dissipation term D. It is 
clear from the Navier-Stokes equations that: 

^ = —diViv"^ - - diVip 
The coarsed-grained expression in the law frequency limit 

S ~ -^^"^"^'^(l + - 4)) (37) 
To arrive in the expression for D we assume Vr ~ Av leading to an expression very similar 
to Kolmogorov's refined similarity hypothesis. 



8 Three-Dimensional Flow 

The most important feature of two-dimensional turbulence, considered in a previous section, 
is irrelevance of the dissipation processes in the inverse cascade range when d < dc- It is this 
irrelevance which was responsible for the gaussian probability density of transverse velocity 
differences. 

The model for the dissipation contribution D in the limit 772 — is evaluated from (37) 
readily. We would like to keep at least some information about Av and the expression must 
be invariant under transformation v — > — v and x — > — x. In addition, the expression must 
be local in physical space. Based on this considerations: 

, . . dAv 
ty pa cidjAuAv— — 
or 

where is a dissipation rate of the "v-component to kitenic energy" Ky — \v^. Locality 
of this model is clear since dj-Av — div{xi) + d2v{x2). The problem is in evaluation of the 
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coefficient c{d) since, in principle, it can be singular at ci = dc- Indeed, it is clear that the 
dissipation term is zero at d < dc- However, the point d — dc separating inverse and direct 
cascade ranges, is a singularity due to infinitly large amplitudes of velocity fluctuations. All 
we can say at this point is that the pdf can be represented as a sum of even and odd functions 
of Av. The symmetric part has a growing with d — dc ^ 0+ width, while the width of odd 
one is 0(1). The behaviour of c{d) in the vicinity of dc is not clear. We feel that it is 0(1) 
ai d > dc and zero at c? < dc- 
Thus, we have: 

D ^ c{d)7]3dr^,dr,,drZ ^ c{d)r]l < AiiA^;-^e^2^"+''3A^ ^ +0{drd^,Z) (38) 

This expression obeys the basic symmetries of the Navier-Stokes equation. The expression, 
similar to (38), was being used in engineering turbulent modelling based on the low-order in 
the Reynolds number coarse-grained equations. 

The last term in the right side of (38), simply modifies the coefficient in front of the first 
term in the left side the equation (6) and does not generate anything new. The expression 
(38) resembles Kolmogorov's refined similarity hypothesis, connecting the dissipation rate, 
averaged over a region of radius r, with (Au)^- Thus, in the limit 772 0: 



[dr,idrj^+-dr,^+{l+h) — -—— — \-c{d)r]3drj^dr,sdr]Z {r]2 = 0,773,r) = —{l-cos{kfr))Z-b^dsZ 
r r Ojj^orj^ 6 rs 

(39) 

The d — dc> counterpart of equation (35) is: 



[d^^ + (l + h + b)^^ + c{d)rjsd,,dr]Z if), = 0, 773, r) = a^^—^^$^rjlZs (40) 

with \E'(?73) chosen in such a way that the generating function Z(0, 0,r) = 1. We consider 

two limiting cases. 

Small — ScaleDynamics. 

Inverse Laplace transform of (40) without right side gives an equation for the pdf P{Av, r): 
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where /3 oc c{d). Since So^s — 0, the coefficients in (41) are chosen to give so,3 = \ Av\^ — a^Pr 
with an undetermined amphtude 03. This is an assumption of the present theory, not based 
on rigorous theoretical considerations. Seeking a solution in a form So,n =< (Av)** >cx r^" 
gives: 

1 + 3/5 1.15 

^„ = — ^ — -n (42 

^ 3(l + /3n) 3(l + 0.05n) ^ ' 

which was derived in [18] together with (3 ~ 0.05. It follows from (42) that: P(0, r) oc r~'^ 
where k — ~ 0.4 for /3 = 0.05. Very often the experimental data are presented as 

P{X,r) where X = V/r'' with 2// = ^ Ri 0.696 for /3 = 0.05. This gives P{X = 0,r) oc 
j,-K+iJ. ^-0.052 compared with the experimental data by Sreenivassan [19]: —n + fi ^ —0.06. 
Let us write P{V, r) = r~'^F{^, r) = r~'^F{Y, r), so that F obeys the following equation: 

Next, changing the variables again —00 < y = Ln{Y) < 00, substituting this into (43) 
and evaluating the Fourier transform of the resulting equation gives: 

BF dF 
(l-P)r— + PK(ik-k'^)F-ikpr— = (44) 
or or 

with the result: F oc r'^('=), where ^{k) — PKf^^^^^Ln{r / L) with r/L « 1. We have to 
evaluate the inverse Fourier transform: 

/oo 
dke-'^ye<^^^ (45) 
-00 

in the limit y — 0(1) and r ^ so that Ln{r/L) —00. The integral can be calculated 
exactly. However, the resulting expression is very involved. Expanding the denominator of 
^{k) gives : 



f 

J —I 



dke-^Hy^'-^)^-''''^^'' (46) 
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and 



with C = V/r^ and Q(r) = A(3Kj^\Ln{r/ L)\. 

To understand the range of vahdity of this expression, let us evaluate < > using the 
expression (47) for the pdf. Simple integration, neglecting 0(/3^) contributions, gives: < 
yn ^a„ Tffif^i^ Q,^^ = (1 + 3P){n — (3{'n? + 2))/3. Comparing this relation with the exact 
result (43) we conclude that the expression for the pdf, calculated above, is valid in the range 
n » 1 and Pn « 1. The properties of the pdf in the range 3 < ^ < 15 are demonstrated on 
Fig. 1 for r/L = 0.1; 0.01; 0.001. The log-normal distribution (48), is valid in a certain (wide 
but limited) range of the V- variation. It is clear from (42) that neglecting the dissipation 
terms {c{d — rfc) oc /9 = 0) leads to ^„ = n/3, i. e. disappearance of anomalous scaling of 
moments of velocity differences. This result agrees with the well-developed phenomenology, 
attributing intermittency to the dissipation rate fluctuations: the stronger the fluctuations, 
the smaller the fraction of the total space they occupy [8], [14]. To the best of our knowledge, 
this is the first work leading to multifractal distribution of velocity differences as a result of 
approximations made directly on the Navicr- Stokes equations. 

To investigate the probability density P{Y, r) in the limit F ^ we introduce an expansion: 

F(Y,r)^Y.CnY'^f2n(r) (48) 

n 

Substituting this into (44) gives: 

f PK2n(2n-l) 

/2nOc(-) (49) 

It is seen from (48)- (49) that the pdf starts bending from the log- normal slope (47) toward 
arF(y,r) = Oaty = Oat: 



r 0K 

Y < (-)^ (50) 
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This inequality shows that as r ^ the pdf develops a narrow cusp at the origin Y = 0. 
If the probability density is plotted in the dimensionless variable X, the bending starts at 
X r° °^. This value was calculated, as above, with P — 0.05. 

Large-scale limit: r/L ^ 1. Now, let us investigate the large-scale limit £ — > 1- Reahzing 
that j3 can be an 0(1) constant we, for the illustration purposes, will investigate the large 
scale limit pretending that j3{d) — >■ as ci — lic- This is also useful since the estimated value 
of /? ~ 0.05 at (i = 3 which is numerically small. In this limit the right side of (41) is 0{rjlZ) 
and cannot be neglected. Repeating the procedure leads to an equation: 

dP ^i±^d d dP _ P d^Z 

dr^ 3r dV^ ^ dV^ dr ~ {Y>r)i dV^ ' 

where a is a proportionality coefficient and r ^ L. As one can see from this equation in the 
limit of small (3 the solution to this equation approaches gaussian. It is also clear that for 
any finite /3, the tails of the pdf are strongly non- gaussian when 

f3Y^ » 1 (52) 

This estimate means that, according to the theory presented above, the perturbative treat- 
ment of deviations from the mean-field gaussian theory is possible but it involves two pa- 
rameters: the ratio e — 1 — j « 1 and P « 1. The fact that the "real-life" P ^ 1/20 may 
explain why the experimentally observed pdf of the large-scale (e << 1 ) velocity fluctuations 
was so close to the gaussian [see [8] and references therein] . 

It is also seen from (52) that at (|^)^ /3 ~ 0.05 the pdf is dominated by a gaussian central 
part. 



9 Conclusions 

The equation (6) formulates theory of turbulence in terms of "only" two unknowns: pressure 

and dissipation terms Ip and D, respectively. It provides a mathematical testing ground for 
various analytic expressions and models obtained from numerical simulations. 
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Armed with the experimental and numerical data, supporting gaussian statistics of trans- 
verse velocity differenced in two-dimensional flows, we showed that the mean field approx- 
imation (the lowest-order term of the expansion (28)) for the pressure contributions (29) 
leads to both Kolmogorov scahng and gaussian statistics of transverse velocity differences. 
In addition, the equation (6) shows that the single-point pdf 's in 2d turbulence are gaussian. 
It is to be stressed that 2d turbulence cannot be a gaussian process and probability density 
P{Au, Av, r) — is not a gaussian. It is only pdf P{Av, r) — J P{Au, Av, r)dAu is a gaussian. 
This statement violates no dynamic constraints. 

One of the most interesting outcomes of the present theory is a discovery of existence of 
the two time-scales in the system which are very different in the vicinity oi d = dc- This 
difference enables one to coarse-grain the Navier-Stokcs equations and neglect all high-order 
non-linearities, generated by the procedure. Using this result the model for for the dissipation 
term D was derived. 

The as yet unresolved ambiguity of this model is its behaviour as d ^ dc- If transition 
from 3d to the non-intermittent state ai d < dc is smooth, then /3 ^ and the resulting 
equation shows onset of both anomalous scaling and non-gaussian statistics. The transition 
can be singular, however: right at d > dc the coefficient P can become 0(1) and a weakly 
intermittent state and weak coupling limit do not exist. In this case, due to existence of 
the small parameter, enabling evaluation of the dissipation expression D, the theory non- 
perturbatively predicts both the shape of the pdf and scaling exponents provided the small 
parameter 5*3^0/ ('S'2,0)^ ^ as d ^ dc- This result is possible since D = Dq + 0{d — dc) 
avd even in the limit d dc, the model D — 0(1). Atd — cic<0,I? = leading to the 
gaussian pdf of transverse velocity differences. Experimental investigation of hydrodynamics 
in a non-integer space dimension is impossible. However, it was demonstarted by Jensen [20] 
that a force-driven shell model is yields the changing sign of the energy flux upon variation 
of a leading parameter. Numerical solution at a critical point (zero flux) demonstrated 
an unsteady state with the growing total energy and the energy spectrum concentrated 
in the vicinity oi kf- The calculat ion also gave kolmogorov energy spectrum at "d > d" 
with growing Kolmogorov constant as d — > dc- It is not yet clear how the eddy viscosity 
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Ln(P(^,r) 




Figure 1: Ln{P{^,r)). Prom bottom to top: r/L=0.1; 0.01; 0.001, respectively. 

approximation works for the shell model, but, since the 0(1) energy flux is flxed by the 
forcing function, the growth of kinetic energy must be related to a relaxation time at 
d = dc and a corresponding small parameter. This result also shows that the phenomenon 
is very robust: all one needs is a point at which energy cascade changes its direction. The 
results of the shell model investigation will be published elsewhere [21]. 
The expression (47) is similar to the one obtained in a ground-breaking paper by Polyakov 
on the scale-invariancc of strong interactions, where the multifractal scaling and the pdf were 
analytically derived for the first time [22]-[23]. In the review paper [23] Polyakov noticed 
that the exact result can be simply reproduced considering a cascade process with a heavy 
stream (particle) transformed into lighter streams at each step of the cascade (fission) . Due 
to the relativistic effects the higher the energy of the particle, the smaller the angle of a 
cone, accessible to the fragments formed as a result of fission. Thus, the larger the number 
of a cascade step, the smaller is the fraction of space occupied by the particles [23] . 
The theory, presented in this paper describes many experimental observations. Still, un- 
derstanding of the limits of vahdity of expression (29) is crucial for the final assessment of 
the theory. The relation (32) shows that (29) is consistent with the gaussian tails of the 
pdf. However, at the present stage we are unable to prove that (29) is the only expression 
leading to this result. The problem is that without experimental detection of at least some 
deviations from the the gaussian statistics of transverse velocity differences, one will not be 
able to understand the limits of vahdity of (29). Given the state-of-the-art of numerical 
simulations, this goal may not be that simple. 
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